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In this paper, we review the recent developments which had taken place in the domain of quasi 
one dimensional Bose-Einstein Condensates (BECs) from the viewpoint of integrability. To start 
with, we consider the dynamics of scalar BECs in a time independent harmonic trap and observe 
that the scattering length can be suitably manipulated either to compress the bright solitons to 
attain peak matter wave density without causing their explosion or to broaden the width of the 
condensates without diluting them. When the harmonic trap frequency becomes time dependent, 
we notice that one can stabilize the condensates in the confining domain while the density of the 
condensates continue to increase in the expulsive region. We also observe that the trap frequency 
and the temporal scattering length can be manoeuvred to generate matter wave interference patterns 
indicating the coherent nature of the atoms in the condensates. We also notice that a small repulsive 
three body interaction when reinforced with attractive binary interaction can extend the region of 
stability of the condensates in the quasi-one dimensional regime. 

On the other hand, the investigation of two component BECs in a time dependent harmonic trap 
suggests that it is possible to switch matter wave energy from one mode to the other confirming the 
fact that vector BECs are long lived compared to scalar BECs. The Feshbach resonance management 
of vector BECs indicates that the two component BECs in a time dependent harmonic trap are more 
stable compared to the condensates in a time independent trap. The introduction of weak (linear) 
time dependent Rabi coupling rapidly compresses the bright solitons which however can be again 
stabilized through Feshbach resonance or by fmetuning the Rabi coupling while the spatial coupling 
of vector BECs introduces a phase difference between the condensates which subsequently can be 
exploited to generate interference pattern in the bright or dark solitons. 

PACS numbers: 


I. INTRODUCTION 

We are aware that eventhough matter pervades the en¬ 
tire universe, it is found in just a few admissible forms 
such as solid, liquid and gas. It is obvious that one can 
initiate a phase transition between different states of mat¬ 
ter by either increasing the temperature or pressure. This 
understanding of generating new states of matter by in¬ 
creasing the temperature was exploited in 1879 by Sir 
William Crookes |1] to create “plasma,” a gas containing 
non negligible number of charge carriers. It must be men¬ 
tioned that the physical states of matter change in going 
from one phase to another while the chemical composi¬ 
tions of matter remains the same. Can one go down the 
temperature scale and generate a new state of matter at 
ultra cold temperatures? This question has been plague- 
ing the minds of scientists and it was in 1995, Cornell and 
Wieman [2] created a “Bose-Einstein condensate” (BEC) 
at super low temperatures. Eventhough envisioned first 
by Albert Einstein and a young Indian physicist named 
Satyendra Nath Bose in the 1920s Ena, it took more 
than seven decades to realize this singular state of mat¬ 
ter. In contrast to plasmas containing superhot and super 
excited atoms, BECs were created at colder and colder 
temperatures near absolute zero and they are composed 
of supercold and super unexcited atoms (see Fig. 0 ). 


* Electronic address: vittal.cnls@gmail.com 



FIG. 1: The energy levels of different physical states of mat¬ 
ter. 


At such low temperatures, a large fraction of the atoms 
get piled up either in the ground state or in the long lived 
metastable state. In other words, the atoms merge to¬ 
gether losing their individual identities and behave like a 
giant matter wave. This phenomenon is known as “Bose- 
Einstein condensation.” 


II. GROSS-PITAEVSKII (GP) EQUATION 

To investigate the dynamics of a Bose-Einstein conden¬ 
sate, we now employ a Hartree or mean field approach 
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and assume that the wave function is a symmetrized 
product of single particle wave functions. In the fully 
condensed state, all bosons (atoms with integral spin) 
are in the same single-particle state, 0(r) and therefore 
one can write down the wave function of the TV-particle 
system as 


N 

^(ri,r 2 ,...,r N ) = n<^). w 

i=l 

The single-particle wave function 0(r i) is normalized 
as 



( 2 ) 


This wave function does not contain the correlations 
produced by the interaction when two atoms are close to 
each other. These effects are taken into account by using 
the effective interaction Uq5(y — r'). According to mean 
field theory, the effective Hamiltonian may be written as 


N 








+ UoY J 8(ri~ r j), ( 3 ) 

i<j 


V{ri) being the external potential. The energy of the 
state given by equation ([l]) is written as 


E = N dr 


/■ 


T—|V^(r)| 2 + F(r)|^(r)| 2 + +- 


: W( r)| 4 

( 4 ) 

From the macroscopic theory for the uniform Bose gas, 
the relative reduction of the number of particles in the 
condensate is of the order of (na 3 ) 1 / 2 , where n is the 
particle density. If we introduce the wave function of the 
condensed state as 


V’(r) = N 1/2 (j)(r), (5) 


then the density of particles is given by 


n(r) = |V’(r )| 2 


(6) 


and for N 1, the energy of the system may therefore 
be written as 


*-/* 


^|V/>(r )| 2 + V>Mr )| 2 + \u 0 m r )| 4 


(7) 

To find the optimal form for t/j, we minimize the energy 
with respect to independent variations of ^(r) and its 
complex conjugate ^(r), subject to the condition that 
the total number of particles 


N = 



(8) 


be constant. For this, one writes SE — ySN = 0 where 
the chemical potential y is the Lagrange multiplier that 


ensures constancy of the particles and the variations of ^ 
and -0* may thus be taken to be arbitrary. Equating to 
zero the variation of E — yN with respect to ^(r) gives 
the following evolution equation m 

- ^~V 2 ^>(r) +V(r)ip(r) + Uo\ip{r)\ 2 ip(j) = mV’W- (9) 

We call equation © as the “time-independent Gross- 
Pitaevskii equation.” This has the form of a time- 
independent Schrodinger equation in which the potential 
acting on particles is the sum of the external potential 
V(r) and a nonlinear term Uo\'i/j(r)\ 2 that takes into ac¬ 
count the mean field produced by the other bosons. 

If one has to look for the dynamics of condensates, it is 
natural to use the time-dependent generalization of the 
Schrodinger equation with the same nonlinear interaction 
term and obtain the time-dependent GP equation 

- Lv^(r) + V(r)if>(r) + U 0 \ip(r)\ 2 ip(r) = ih 

( 10 ) 

To ensure consistency between the time-dependent GP 
equation ( [To| ) and the time-independent GP equation ©> 
under stationary conditions ^(r, t) mu st evolve in time 
as exp(— iyt/K). 

In the above equation ( |l0| ), ip(r , t), r = (x,y,z) rep¬ 
resents the condensate wave function, V 2 denotes the 
.Laplacian operator, V(r) is the trapping potential as¬ 
sumed to be V(r) = m(u; 2 r 2 +cj 2 x 2 ) where r 2 = y 2 + z 2 , 
uj r , x are the confinement frequencies in the radial and ax¬ 
ial directions respectively, Uo = 47r h 2 a/m corresponds to 
the strength of interatomic interaction between the atoms 
characterized by the short-range 5-wave scattering length 
a, and m is the atom mass. 

From equation ( fIo| ) it is obvious that the GP equa¬ 
tion is an inhomogeneous (3+1) dimensional nonlinear 
Schrodinger (NLS) equation. The inhomogeneity origi¬ 
nates from the potential V(r) that traps the atoms in 
the ground state and the nonlinearity coefficient Uo which 
represents the interatomic interaction related to the scat¬ 
tering length a. The scattering length can have either 
positive or negative values which means the interaction 
can be either repulsive or attractive. It should be men¬ 
tioned that it should be possible to vary the scattering 
length periodically with time a(t) employing Feshbach 
resonance [8]. This means that understanding the dy¬ 
namics of BECs boils down to solving a variable coeffi¬ 
cient (3+1) NLS equation for suitable choices of trapping 
potentials V(r) and temporal scattering lengths a(t). 


III. INTEGRABILITY OF GP EQUATION- 
REVIEW OF ANALYTICAL METHODS 


Looking back at the (3+1) dimensional Gross- 
Pitaevskii equation (10), it is obvious that it is in gen¬ 
eral nonintegrable for an arbitrary trapping potential and 
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interatomic interaction. Hence, one has to investigate 
whether the GP equation would admit integrability in 
lower spatial dimensions for specific choices of trapping 
potentials and scattering lengths. In other words, one has 
to look for the associated nonlinear excitations in (1+1) 
and (2+1) dimensional GP equations which would reflect 
upon the integrability of the associated dynamical sys¬ 
tem under consideration. In a three dimensional BEC, 
when the transverse trapping frequency uo r (r = y , z) 
is very high compared to the longitudinal trapping fre¬ 
quency uj x , then the transverse confinement is too tight to 
allow scattering of atoms to the excited states of the har¬ 
monic trap in the transverse direction. Under this con¬ 
dition, one obtains “cigar” shaped BECs and the three 
dimensional GP equation becomes quasi one-dimensional 
in nature. Again, it should be mentioned that the quasi 
one-dimensional GP equation can be shown to be inte¬ 
grate only for suitable choices of trapping potential V ( x ) 
and interatomic interaction Uq(x). 

A. Analytical Methods 

Eventhough one cannot precisely define the concept of 
integrability of a dynamical system governed by a nonlin¬ 
ear partial differential equation (PDE), one can look for 
the possible signatures of integrability namely, Painleve 
(P-)property [9;, Lax-pair m soliton solutions etc. A 
given nonlinear dynamical system governed by a nonlin¬ 
ear PDE is said to admit P- property if the corresponding 
solution can be locally given in terms of a Laurent se¬ 
ries expansion in the neighborhood of a movable singular 
point/manifold. The existence of a Lax-pair of a given 
nonlinear pde implies that one can somehow linearize the 
nonlinear dynamical system and subsequently exploit it 
to generate soliton solutions, thereby consolidating its in¬ 
tegrability. In this section, we dwell upon the analytical 
techniques like Inverse Scattering transform PH. Gauge 
transformation method m, Darboux transformation ap¬ 
proach [IMS], Hirota’s direct method [16], S3 besides 
the approximation method like variational approach. 

I. Inverse Scattering Transform 


The Inverse scattering Transform is a nonlinear 
analogue of the Fourier transform which has been em¬ 
ployed to solve several linear partial differential equa¬ 
tions. Given the initial value of the potential q(x, 0) and 
the boundary conditions, one has to identify two linear 
differential operators L and B so that one can convert 
a (1+1) dimensional nonlinear pde into two linear equa¬ 
tions, namely a linear eigenvalue problem 

L<f> = A<f>, (11) 

and a linear time evolution equation 

&t = ( 12 ) 


such that the compatability condition of the above two 
equations © and ( [T2] ), z.e., L t = [B,L\ generates the 



FIG. 2: Schematic diagram of the inverse scattering transform 
method. 


nonlinear PDE one has started with. Once the lineariza¬ 
tion is performed in the above sense for a given nonlinear 
dispersive system q t = K(q), where K(q) is a nonlinear 
functional of q and its spatial derivatives, the Cauchy ini¬ 
tial value problem corresponding to the boundary condi¬ 
tion q —)> 0 as \x\ -+ 00 can be solved by a three step 
process indicated schematically in Fig. [2] 

This method involves the following three steps: 


1. Direct scattering transform analysis 

Considering the initial condition q(x, 0) as the po¬ 
tential, an analysis of the linear eigenvalue problem 
0 is carried out to obtain the scattering data 
5[ 0). For example, for the KdV equation we have 

5(0) = A n (0), n = 1 , 2 ,..., TV, C n (0), R(x , 0), —00 < x < 00 

(13) 

where N is the number of bound states with eigen¬ 
values A n , C n ( 0) is the normalization constant of 
the bound state eigenfunctions and R(x, 0) is the 
reflection coefficient for the scattering data. 


2. Time evolution of scattering data 

Using the asymptotic form of the time evolution 
equation (12) for the eigenfunctions, the time evo¬ 
lution of the scattering data S(t) can be deter¬ 
mined. 


3. Inverse scattering transform (1ST) analysis 

The set of Gelf and — Levitan — Marchenko in¬ 
tegral equations corresponding to the scattering 
data S(t) is constructed and solved. The resulting 
solution consists typically of N number of local¬ 
ized, exponentially decaying solutions asymptoti¬ 
cally (t ± 00 ). In this way, one can success¬ 
fully solve the initial value problem of the nonlinear 
PDE. 


From the above, it is obvious that solving the initial 
value problem of the given nonlinear PDE boils down to 
solving an integral equation. 
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II. Gauge Transformation Approach 


This is an iterative method enabling one to gen¬ 
erate soliton solutions starting from a seed solution. 
In this method, one again begins with the Zakharov- 
Shabat(ZS)-Ablowitz-Kaup-Newell-Segur (AKNS) linear 
systems PUCE] given by the following equations: 

= U$ (14) 

= V§ (15) 

where 

U =(~r X l I) • V ~(cD )■ <“> 

In the above equation, that is equation ( p~6] ) , A repre¬ 

sents the spectral parameter, (g, r) the potential/field 
variable of the given nonlinear PDE, A, B , C and D cor¬ 
respond to the undetermined functions of A, x, £, g, r and 
their spatial/time derivatives. The compatibility condi¬ 
tion ($ x ) t = ($ t ) x yields U t - V x + [U,V] = 0 which 
is equivalent to the given nonlinear PDE. One generally 
calls the matrices U and V as “Lax-pair” which contains 
information about the linearization of the given dynam¬ 
ical system. It should be mentioned that for scalar non¬ 
linear PDEs, U and V are 2x2 matrices and the eigen¬ 
function <f> is 2x1 column vector while for vector (two 
component in general) nonlinear PDEs, the eigenfunc¬ 
tion <f> is a 3x1 column vector and U and V are 3x3 
matrices. 

Now, gauge transforming the eigen function <f>, we have 
an iterated eigen function &A) = g <£ where g = g(x,t, A) 
is a 2x2 matrix function so that 

$(i) = j7(i)$ (17) 

(is) 

where, 

u w = gUg- 1 + g x g-\ (19) 

W =gVg- 1 +g t g- 1 . (20) 

The transformation function g must be adopted from the 
solutions of certain Riemann problems in the complex A 
plane and it must be a meromorphic (regular) function. 
The simplest form of g satisfying the above criteria can 
be written as 


9 = 


/ + 


Ai — hi 

A -Ai 


P(x,t) 



( 21 ) 


where Ai and Hi are two arbitrary complex numbers and 
P is an undetermined 2x2 projection matrix (P 2 = P). 
Hence, g~ x is now given by 


9 


-l 



0 

-1 



Ai — Hi 
X — Hi 


P{x,t) 


( 22 ) 


Thus, the vanishing of the apparent residues at A = Ai 
and X = Hi imposes the following constraints on P as 


Px = (/ - P)cr 3 P(/ii)cr 3 P - Pcr 3 P(A 1 )cr 3 (/ - P)(23) 
Pt = (/-P)a 3 y(/ii)a 3 P-Pa3E(A 1 )a3(/-P)(,24) 

where 


^=(j_°i)- (25) 

To generate a new solution from a given solution (vac¬ 
uum), for example g° and r° associated with matrices U° 
and V° , the eigen value problem takes the following form 

$0 = $0 = V 0$0' (26) 

where U\ see d = U° and V\ see d = U°. Now, one can solve 
equations ( [23] ) and ( [24] ) using the vacuum eigen function 
4>° such that 


P = cr 3 Pcr 3 , (27) 


where 


P = / [trace 

and is a 2x2 matrix defined by 


(28) 


i) ( Ai)" 1 , (29) 


where m i and n\ are arbitrary complex constants and 
4>° is a solution of the vacuum linear system governed by 
equations (26). 

Now, substituting the meromorphic function g with the 
projection matrix P given by equation (27) into equation 
(19), we get 


V'W = (:$ '/)-*<*■-»)( X P d 2 )- 

(30) 

A comparison of the eigenvalue problems expressed 
in terms of the vacuum eigenfunction <f>° and the new 
(transformed) eigenfunction 4>( 1 ) would enable us to re¬ 
late the vacuum solution of the associated nonlinear PDE 
with the new solution. Thus, we get an explicit solution 
as [12] 


q W = - q o _ 2*(A! - A*i)Pi2, (31) 

r (1) =-r° + 2i(Ai-A4i)P2i, (32) 


where 


Pl2 = 


M\ 


(i) 


12 


M 


(i) 


li 


■ M< 


(i) : 


22 


P21 = 


Mr 


( 1 ) 


21 


M 


(i) 


li 


Mr 


( 1 )‘ 


22 


(33) 


Once we get q W and A 1 ^ from the input solution g° 
and r°, one can repeat the same procedure to obtain yet 
another new solution q^ and A 2 ^) using q W and A 1 ^) 
as the input solution. For example, to construct second 
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iterated solution q^ 2 \ r^ 2 \ one needs to find a solution 
of the linear system given by equations ( [17} and (18) 
where the matrices U 1 and V 1 are associated with the 
input solutions gD) and r^ 1 ). However, one can find the 
solution of above eigen value problem in terms of <L° as 


Now, introducing the Darboux transformation into the 
known eigen function, we now obtain the transformed 
eigenfunction as 

$ (1) (x, t, A) = D(x, t, A)$(z, t, A), (43) 


= g&°. 


(34) 


Thus, the new iterated solution q^ and r^ (in anal¬ 
ogy with equations (|3Tj) and (32)) can be written as 


q {2) = -q {1) - 2«(a 2 - /i 2 )A ( 2 ) . ( 35 ) 

r ( 2 ) = _ r (i) _|_ 2i(A 2 - (36) 


Thus, one can repeat the same procedure 7V-times to 
obtain the N th iterated solution as 


q N =-q( N ^ -2i(\ N - /i n )P ^2 1} , (37) 

r N = -r (JV_1) - 2 z(Ajv - (38) 

The above iteration could become extremely handy, 
particularly in the context of the generation of multisoli- 
ton solutions as one can obtain TV-soliton solution from 
the vacuum eigen function 4>° of the linear system. 

III. Darboux Transformation Approach 


where D(x, £, A) is the “Darboux matrix” which is equiv¬ 
alent to XI—S. While I is the identity matrix, the matrix 
S can be generated as 

S = HAH ~\ detH/0, (44) 


where the matrix H is defined as H — (hi,..., hjy), 
where hi represents the column solution of the linear 
eigenvalue problem given by equations ( fl4| ) and ( p~5] ) , 
A = diag(Ai,..., An)- The new eigen function 4 A 1 ) again 
satisfies equations (fTTl and (18) and the Darboux matrix 
D plays the role of the transformation function g. Ac¬ 
cordingly, we have 


[/W = DUD- 1 + D x D~ l , (45) 

V {1) = DUD- 1 + D t D~ l . (46) 


To check the form of U' l> given by equation (45), we 
substitute eqns. into the x derivative of equa¬ 

tion^) to obtain 


In 1882 G. Darboux studied the eigen value problem of 
the one dimensional Schrodinger equation 


- ®xx ~ = A4>, 


(39) 


where q{x) is a potential function and A is a constant 
spectral parameter. He postulated that if q(x) and 
4>(x,A) are two functions satisfying equation (39) and 
f(x) = 4>(x,Ao) is a solution of equation of (39) for 
A = Ao where Ao is a fixed constant, the functions q' 
and 4>' are defined by 


q' = <? + 2 (ln/) a 


(40) 


with 


& xx - q'& = a $'. 


(41) 


From equations (39) and (Dll), it is obvious that they 


are of the same form. Therefore, the transformation (40) 


which converts the functions (g, <f>) to (^ / ,4> / ) satisfying 
the same partial differential equation is the original Dar¬ 
boux transformation which is valid for / ^ 0 . 

In this method, one again begins with the ZS-AKNS 
linear eigen value problem given by equations (14) and 


(15), where the eigen function is a 2 x 2 matrix of the 


following form 

V ' $22OM, A) 


(42) 


t/ (1) $ (1) = D X $ + DU$. 


(47) 


Then, making use of equation (43), one obtains equa¬ 
tion (45). Similarly, one can check the form of given 
by (46}. 

Thus, starting from a seed solution U of the given non¬ 
linear pde, one can generate the iterated U^\ This pro¬ 
cedure can be repeated to generate multisoliton solutions. 


IV. Hirota Bilinear Method 


Although the inverse scattering formalism was the first 
analytical technique that has been developed to solve 
the initial value problem of nonlinear pdes, it involves 
sophisticated mathematical concepts like solving an inte¬ 
gral equation and hence quite complicated and intricate. 
Moreover in this method, one should have a prior knowl¬ 
edge of the potential u(x,t) at t = 0 , namely the ini¬ 
tial data u(x, 0 ) and the boundary conditions imposed 
on it. On the other hand, eventhough Darboux and 
gauge transformation approaches are iterative in nature 
and are purely algebraic without involving highly com¬ 
plex mathematics, they warrant the identification of the 
Lax-pair of the associated dynamical system. Hence, it 
is imperative to look for an alternative method to gen¬ 
erate localized solutions (soliton solutions) of nonlinear 
pdes and in this context, Hirota’s direct method comes 
to our rescue. In this method, neither does one need any 
prior information about the potential (or physical field) 
of the associated nonlinear pde, nor the Lax-pair of the 
associated dynamical system. This method which has an 
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inbuilt deep algebraic and geometric structure is more el¬ 
egant and straightforward and can be directly employed 
to generate soliton solutions of nonlinear pdes. 

The salient features of the Hirota method are the 
following: 

i) The given nonlinear partial differential equation 
has to be converted into a bilinear equation through a 
transformation which can be identified from the Painleve 
analysis. Each term of the bilinear equation has the 
degree two. 

The Hirota bilinear operators are defined as 


alizable potentials. To start with, we consider the dy¬ 
namics of BECs in a time independent harmonic trap 
with exponentially varying scattering length. We gener¬ 
ate the associated bright solitons and study their colli- 
sional dynamics. We then introduce suitable time depen¬ 
dence in the harmonic trap and investigate the dynamics 
of BECs. We then show that how the interplay between 
trap frequency and temporal scattering length can gen¬ 
erate matter wave interference pattern in the collision of 
bright solitons. We then reinforce the binary attraction 
with a repulsive three body interaction to enhance the 
stability of BECs. 


D?D n x {G-F) 


f d_ _ d_\ m ( d_ _ d_\ n 

dt') \3x dx') 

G(t,x)F(t',x')\t>= t ,x'=x 


ii) The dependent variables G and F in the bilinear form 
have to be expanded in the form of a power series in 
terms of a small parameter 6 as 


G = ££ (1) + eV 3) + e:V 5) H-, (48) 

F = 1 + e 2 / (2) + e 4 / (4) + ''' • (49) 


A. Dynamics of quasi one dimensional BECs in a 
time independent harmonic trap 

For a time independent harmonic oscillator potential 
and exponentially varying scattering lengths, the GP 
equation ( [To| for cigar shaped BECs takes the following 
form p20l424j 

+ 2 a(t) \ip \ 2 V> + ^A 2 x 2 ip = 0, (50) 


iii) After substituting the above functions into the bilin¬ 
ear form and equating different powers of £, a set of linear 
pdes can be generated. 

iv) Finally, solving the linear pdes, one can generate soli¬ 
ton solutions. 

It should be mentioned that the key to the success of 
the method lies in the identification of the dependent 
variable transformation as well as in choosing an opti¬ 
mum power series to linearise the given nonlinear pde. 


A. Approximation method 
I. Variational Approach 

Variational approach is a qualitative semi-analytical ap¬ 
proach. By using the variational approximation, one 
can study the dynamics of Bose-Einstein condensates 
described by the mean-field Gross-Pitaevskii equation. 
For the purpouse of variational analysis, Larangian den¬ 
sity is calculated for the corresponding time-dependent 
Gross-Pitaevskii equation and the effective Lagrangian 
can be obtained by integrating the initial trial wave 
function with variational parameters over space. The 
numerical value of each variational parameter can be 
obtained from the numerical solution of corresponding 
Euler-Lagrangian equations. 


IV. BRIGHT MATTER WAVE SOLITONS AND 
THEIR COLLISION IN SCALAR BECS IN 
CERTAIN SIMPLE POTENTIALS 


where time t and coordinate x are measured in units 
2/w± and a_ j_, where a± = ( h/muj^) 1 / 2 and ao = 
(h/muj o) 1 / 2 are linear oscillator lengths in the transverse 
and cigar-axis direction respectively. uo±_ and ujq are 
corresponding harmonic oscillator frequencies, m is the 
atomic mass and the trap frequency A = 2\ujq\/uj±_ << 
1. The Feshbach resonance managed nonlinear coeffi¬ 
cient which represents the scattering length reads a{t) = 
doexp(At). 


Equation (50) admits the following Lax-pair 


V 


= £/$, U = 


i( Q 
-Q* -< 


(51a) 

V$, (51b) 

-2 i( 2 + i\x( + i \Q\ 2 [(Ax - 2 QQ + iQ x \ 
[-(Ax - 2()Q* + iQ*\ 2i( 2 — i\x( — i \Q\ 2 


where we have slightly modified the Lax-pair (given in 
ref. 23 ) by allowing the nonisospectral parameter ( to 
be complex keeping the initial scattering length unity. 
The nonisospectral complex parameter obeys the first or¬ 
der ordinary differential equation of the form 


Ct = AC, CW = ^) + W, (52) 


and the macroscopic wave function ^ is related to Q by 
the transformation 


Q = exp 




(53) 


In this section, we investigate the dynamics of scalar 
Bose-Einstein condensates in certain simple physically re- 


Now, to generate the bright soliton solutions of equa¬ 
tion (50), we consider a trivial vacuum solution =0 
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to give the following vacuum linear systems 


$(°) = ^ ^ ^ $ (0) = P (0) $ (0) , (54a) 

$(o) = (-Zi^+iXxC 0 \ $( o) 

1 y 0 2 iC? — i\x( J 

= V^^ 0) . (54b) 


Solving the above linear systems keeping in mind that 
the spectral parameter ( varies with time by virtue of 
equation ( [52] ) ? we have 


$ (o) 0 x,t,0 = 


g ixC,—2i f C dt Q 

q e -ixC+2i J C 2 dt 


(55) 


Now, effecting the gauge transformation 

$ ( 1 ) (M,C) ~s$ ( 0 ) OM,C), ( 56 ) 


where “g” is a meromorphic solution of the associated 
Riemann problem. The new linear eigenvalue problems 
now take the following form 

$<4) = = yW^ 1 ), (57) 

with 

u {l) = gU^g- 1 + g x g~ l , (58a) 

=gV^g- 1 +g t g- 1 . (58b) 

We now choose g as 

9 = ( x + ^a Pl{x,t) ) (o -1 ) (59) 

where Cl and fii are arbitrary complex parameters and 
Pi is a projection matrix (P 2 = Pi). Imposing the con¬ 
straint that U W and do not develop singularities 
around the poles C = Ci and ( = H 1 , the choice of the 
projection matrix Pi is governed by the solution of the 
following set of partial differential equations 

P lx = (l-P 1 )a 3 U^(g 1 )a 3 P 1 -P 1 a 3 U^ 0 \Ci) 

ct 3 (1 - Pi) (60a) 

P u = (1 - - P 1 a 3 V^°\Ci) 

(7 3 (1 - Pi) (60b) 

where 


03 ( 0 -1 ) ' ( 61 ) 

Looking at the above system of equations, we under¬ 
stand that Pi depends only on the trivial matrix eigen¬ 
function 4>( 0 )(x, £, £), a diagonal matrix and has a com¬ 
pact form given by 


Pi(x,t) 


^3 


77-VTTiyr^s, (62a) 

[trace MvBJ 

n! (62b) 


L 

FIG. 3: The dynamics of bright solitons for the parametric 
choice A= 0.02, /3 0 =2, o 0 =0.1, <h=0.5, 0i=O.l. 



FIG. 4: Contour Plot of bright soliton for A=0.02. 


where mi and ni are arbitrary complex constants. 
Hence, choosing the complex parameters (j = aq (t) + 
i/3i(t) and = Q and employing the gauge transforma¬ 
tion approach [ 12 ], we arrive at the matter wave bright 
soliton 


ip W (x,t) = 2/3 0 exp sech(8i)exp(if;i), 

(63) 

where 

6 >i = 2fi\x — 8 J(a 1 p 1 )dt + 26 1 , (64a) 

Ci = 2a\x — 4 J (af — Pf)dt — (64b) 

«i = ai 0 e At , Pi = /3 w e xt . (64c) 


The ab ove brigh t soliton solution given by equations 
(63) and ( 64af 64c ) is identical to the one given by Liang 
et al. m using Darboux transformation. From the 
profile of the bright soliton trains shown in Fig. [3] we 
infer that matter wave density |^| 2 increases with the 
increase of the absolute value of the scattering length 
leading to the compression of the bright soliton trains of 
BEC. The contour plot (shown in Fig. [4| which takes a 
cross-sectional view of Fig. [3] in the x — t plane shows 
that the width of bright solitons decreases progressively. 



FIG. 5: The dynamics of bright soliton for A= -0.02. 









2DC 



FIG. 6: Contour plot of bright soliton for A= -0.02. 


Figure [5] shows that the peak value of the matter wave 
density |^| 2 decreases with the decrease of the absolute 
value of the scattering length leading to a broadening of 
the bright soliton trains thereby enhancing their width 
and this is again confirmed by the corresponding contour 
plot in figure [ 6 ] Our investigation shows that the scatter¬ 
ing length can be suitably manipulated to compress the 
bright solitons of BECs into an assumed peak matter den¬ 
sity without causing their explosion while on the other 
hand, it can be manoeuvred judiciously to broaden the 
localized solitons without allowing the dilution of the con¬ 
densates and this interpretation completely agrees with 
that of Liang et al. [ 22 ]. Investigations of the quasi one 
dimensional GP equation [ 22 | [23] 1251129] in the presence 
of an expulsive parabolic potential for positive scattering 
lengths also confirm our above observations. 

From the above, we observe that one can either com¬ 
press or broaden the bright solitons in the expulsive time 
independent trap either by exponentially increasing or 
decreasing the scattering lengths respectively. It must 
be mentioned that the exponentially varying scattering 
length with a trap frequency dependence would make the 
exact solutions of the GP equation less interesting from 
an experimental point of view. Hence, it would be in¬ 
teresting to investigate the impact of a general time de¬ 
pendent scattering length and a time dependent trap on 
the condensates. The addition of time dependent trap 
frequency will facilitate us to tune the trap suitably and 
study its impact on the condensates. 

B. Impact of transient trap on BECs 


The introduction of time dependance in the trap en¬ 
sures that the condensates are now confronted with both 
time dependent scattering length and time dependent 
trapping potential and accordingly equation (50) gets 
modified as (in dimensionless units) [30] 


.8^ ld 2 ^ /m,| 2, A(t) 0/ ~ ,nr\ 

+ 2d^ + ^ = °» ( 65 ) 

where a(t) = —2 a s (t)/aB, A (t) = cjq( t)/tJj_, a# is the 
Bohr radius, A(£) describes time dependent harmonic 
trap which can be confining (A(t) < 0) or expulsive 
(A (t) > 0). This equation has been mapped onto a lin¬ 
ear Schrodinger eigenvalue problem m and one soliton 


solution has been expressed in terms of doubly periodic 
Jacobian elliptic functions. However, the above approach 
cannot be employed to generate multisoliton solutions 
analytically. The gauge transformation approach comes 
in handy at this juncture as it offers the advantage of 
constructing multisoliton solution from the solution of 
the corresponding vacuum linear system. 

Now, to generate the bright solitons of equation (65) 
for both regular and expulsive potentials, we introduce 
the following modified lens transformation [29H32] 


^(x,t) = yj A(t)Q(x , £)exp(i<F(x, t)), ( 66 ) 

where the phase has the following simple quadratic form 

§(x,t) =-\c(t)x 2 . (67) 


Substituting the modified lens transformation given by 
equation ( 66 ) in equation (65), we obtain the modified 
NLS equation 


iQt 

with 

and 


\Qxx - ic(t)xQ x - ic(t)Q + a(t)A(t)\Q\ 2 Q = 0 , 

( 68 ) 

(69) 


\(t) =c'(t) -c(t) 2 , 


c(t) = - —lnA(t). 


(70) 


Equation ( 68 ) admits the following linear eigenvalue 
problem 


V = 


where 


= u<l>, 

T T _ ( K{t) 

u - V -Q* 

Q \ 

-m) ’ 

(71) 

= Vct>, 



(72) 

f V n 

v 12 \ 



V^21 

v 22 J 



: -m 2 

+ ic(t)xQ(t ) + 

l -a{t)A(t)\Q \ 2 


■ ( c(t)x - 

C (t))Q + 2 Qx 



■ —(c(t)x 

-mw + \ 

Ql 


: m 2 - 

ic(t)x((t) - l - 

a(t)A(t)\Q \ 2 

(73) 


In the above linear eigenvalue problem, the spectral 
parameter which is complex is nonisospectral obeying 
the following equation 


c\t) = cmt), 


(74) 
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with a(t) = 1/A(t). It is obvious that the compatability 
condition (<t> x )t = (0t) x generates equation ( [68]) . 

Substituting equation (70) with a(t) = 1/A(t) in equa¬ 
tion (69), we get 


a" (t)a(t) - 2a\t) 2 - \{t)a{t) 2 = 0. (75) 

From the above, it is evident that the GP equation ( [65] ) 
is completely integrable only if the trap frequency A Jt) 
and the scattering length a(t) are connected by equation 
(75) and the above condition is consistent with Ref. [29] . 


Thus, the modified lens transformation has facilitated the 
identification of integrability of equation (65). It should 


also be mentioned that equation ( [65] ) is completely inte¬ 
grable only for certain suitable choices of trap frequency 
A (t) depending on the solvability of equation (69). For 
example, when A (t) = constant = k and a(t)=e^q where 
A is the trap frequency, equation ([65]) reduces to (50) 


describing the dynamics of BECs moving in an expulsive 
parabolic potential and exponentially varying scattering 
length [2J] 24 . The above parametric choice is consis¬ 
tent with equation ( [75] ) with k = — A 2 ensuring the inte¬ 
grability of the model. The above model has also been 
experimentally realized [33] . 

To generate the soliton solution of equation ( |68| ) (or 
equation (65)), we consider the seed solution = 0 
and solve the linear systems given by equations <03 and 
(72) keeping in mind equation (74) to obtain 


(76) 

Employing the gauge transformation approach and 
choosing (0 = aq (t) + i0i(£) and /r i = Ci, one obtains 
the one soliton solution of equation (65) [34j 

V(x,t)= .LL-2/h(t) sech#i e -5 c W x2+i « 1 , (77) 

V a W 


where 


0± = 2/3i(t)x-4 f (ai(t)/3i(t))dt + 2£i, 

Jo 

2ai(t)x — 2 f (ai(t ) 2 — 0i(£ ) 2 )dt — 20i(78) 
Jo 


Cl = 


ol\ = a. x 0 efi c(t ' )dt ', 

Pi = Aoe/o^')*', 

and 0i, S i, aio and 0io are arbitrary real constants. 

The striking feature of this bright soliton solution is 
that its amplitude relies strongly on the scattering length 
a(t) and the time dependent trap A (t) while the velocity 
is governed by the external trap A (t) alone. 

This procedure can be easily extended to generate mul- 
tisoliton solution and one can study the collisional dy¬ 
namics of bright solitons for a suitable choice of A (t) and 
a(t) consistent with equation (75). 



FIG. 7: Two soliton interaction in the expulsive trap (A (t) < 
0) with a(t) =aoexp(— 0.125t 2 ), ao — 0.5, aio = 2.31, 0io = 
1.5, 020 = 3.12, 020 = 1.2, 0i = .005, Si = 0.002, 0 2 = 0.002, 
S 2 = 0 . 001 . 



FIG. 8: Two soliton interaction in the expulsive trap (A (t) < 
0) with a(t) =aoexp(— 0.125t 2 ), ao = 0.5, aio = 2.31, 0io = 
1,5, 020 = -2.12, 020 = 1.2, 0i = .05, S 1 = 0.02, 0 2 « 0.02, 
62 = 0.01. 


Thus, it is obvious that one can obtain varieties of 
soliton profiles depending on the choice of the scattering 
length a(t) and the time dependent trap A (t) consistent 
with equation ( [75] ) . 

Figures 0 and 0 describe the evolution of the two 
soliton solution for an expulsive trap (A (t) < 0) for dif¬ 
ferent initial conditions evolving the scattering length of 
the form a(t) = 0.5 exp(—0.125t 2 ). From the figures, one 
observes that the matter wave density |0| 2 of the con¬ 
densates decreases slowly by virtue of the decrease in the 
absolute value of the scattering length and the trajec¬ 
tory of the soliton pulses is dictated by the initial con¬ 
ditions. It should be mentioned that the identification 
of this critical parametric regime in which one observes 
the slow decay of the condensates enables one to avoid 
this domain by operating the system under a safe range 
of parameters. 

Figure[9](a) shows the interaction of solitons for a(t) = 
0.5exp(0.0025£ 2 ). It can be observed from fig. that 
the confining nature of the trap (A (t) > 0) is preserved 
only for a finite length of time (t < 14). During this 
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(a) 



FIG. 9: Two soliton interaction in the confining trap (A (t) > 
0) at different intervals of time with a(t) =ao exp(0.0025t 2 ), 
a 0 = 0.5, aio = 0.01, (3 10 ~ 0.1, a 20 = 0.28, f3 20 = 0.11, 
(f)i — S 2 — 0.1, (f )2 — — 0.2. 


period, the two soliton pulses slide over each other like 
liquid balls as shown in figures [9] (c) and [9] (d). After 
this critical period (t >14), the trap becomes expul¬ 
sive again which sets in the compression of the soliton 
pulses resulting in the increase of the matter wave den¬ 
sity |^| 2 of the condensates. It can also be observed that 
for a(t) = 0.5 exp(—0.125£ 2 ), the absolute value of the 
scattering length decreases and hence there is a slow de¬ 
cay of the condensates while for a(t) = 0.5 exp(0.0025t 2 ), 
eventhough the absolute value of the scattering length 
increases, the soliton pulses begin to get compressed (or 
the matter wave density |^| 2 increases) after a finite time 
delay. It can be easily understood that this delay is in¬ 
troduced by the time dependent trap. When the time de¬ 
pendent trap A (t) becomes a constant, equation (65) re¬ 
duces to the dynamics of BECs in an expulsive parabolic 
potential and time independent scattering length. Un¬ 
der this condition, the soliton trains begin getting com¬ 
pressed and the matter wave density |^| 2 increases as 
soon as the absolute value of the scattering length in¬ 
creases and one does not observe any time delay in the 
compression of soliton trains [221 HU • It should also be 
mentioned that this time delay in the compression of the 
soliton pulses can be suitably manipulated by changing 
the trapping coefficient A (t). Thus, the bright solitons 
can be compressed into a desired width and amplitude in 
a controlled manner by suitably changing the trap and 
our observation is consistent with the numerical results 



FIG. 10: Two soliton interaction in the confining trap (A (t) = 
0.09) with a(t) =ao exp(0.3zt), ao = 0.5, a±o = 0.09, /3io = 
0.71, 020 - 0.031, £20 = 0.11, = 5.1, = 7.2, 0 2 = 4.2, 

S 2 = 4.1. 


m- 

It can also be observed that at t >14, the time depen¬ 
dent trap becomes expulsive again for the same choice 
of a(t) (i.e., a(t) = 0.5exp(0.0025t 2 )) (fig. [9J1). In or¬ 
der to sustain the confining nature of the trap (A (t) >0), 
the scattering length should become complex as it is be¬ 
ing done in the case of cold alkaline earth metal atoms 
[36] . Under this condition, the matter wave density l^l 2 
periodically changes with time by virtue of the periodic 
modulation of scattering length [37l [38 ] and this is remi¬ 
niscent of the recent experimental observation of Faraday 
waves [32]. We also observe that the two soliton pulses 
keep exchanging energy among themselves continuously 
during propagation as shown in Fig. [lO] 

We observe from the above that the addition of time 
dependence in the trap enables one to stabilize the con¬ 
densates for a longer period of time by selectively tun¬ 
ing the trapping potential. It should be mentioned that 
one can also selectively choose A (t) and a(t) and observe 
their interplay in the collisional dynamics of bright soli¬ 
tons. The interplay between A (t) and a(t) consistent with 
equation ([75]) results in the “matter-wave interference 
pattern” in the collisional dynamics of bright solitons. 


C. Matter wave interference pattern in the 
collision of bright solitons 

To generate the matter wave interference pattern, we 
now allow the two bright solitons to collide with each 
other in the presence of a trap for suitable choices of 
scattering length a(t) and trap frequency A (t) (or c(t)) 
consistent with equation ( [75] ). 

Case (i): When c(t) = —1, the scattering length 
evolves as a(t) = aoe _t (shown in Fig. 12 [b)) where ao 
is an arbitrary real constant and the trap frequency A (t) 
becomes expulsive and is equal to a constant (A (t) = — 1). 
Under this condition, the collisional dynamics of two 
bright solitons which are initially separated as shown in 

(a) in the expulsive harmonic trap is shown in 

(b) and the corresponding density evolution in 


Fig. 

Fig. 
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(b) 



FIG. 11: (a) Initial position of two bright solitons at t = — 3 
for c(t) sb —1. (b) Interaction of two bright solitons form¬ 
ing the interference pattern for the choice a(t) = aoe~ f with 
a lo =0.01, o 20 =0.8, #i O =0.06, # 2O =0.012, a o =0.2, 0i =0.01, 
02=0.1, (Ji=0.1, <* 2 =0.01. 




-60 -40 -20 0 20 40 60 



FIG. 12: (a)Contour plot of figJTTJb) depicting matter wave 
density evolution, (b) The time evolution of scattering length 
a(t) corresponding to case (i). 


Fig. [l2] (a). The density evolution consists of alternating 
bright and dark fringes of high and low density respec¬ 
tively while the phase difference between the condensates 
continuously changes with time as shown in Fig. [l4| (a). 
Case(ii): When A (t) = 0.3 — 0.09£ 2 (shown in Fig. 
a)) and a(t) = aoe^ 0-15t \ the soliton interaction and 
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the corresponding contour plot showing the interference 
pattern in the confining region are shown in Figs. [I3| (b) 
and [T3 ] (c) respectively. From the Figs. [l3| (a)-(c)7it is 
clear that when the trap A (t) enters the confining region 
from the expulsive domain, the interaction of the soli¬ 
tons generates the interference pattern in the confining 
regime and the pattern disappears once the trap becomes 
expulsive again. The above choice of A (t) and a(t) can be 
synthesized easily under suitable laboratory conditions. 
The phase difference between the condensates continu¬ 
ously changes with time as shown in Fig. [14] (b) 


Thus, our results reinforce the fact that the matter 
waves originating from the condensates (bright solitons) 
do interfere and produce a fringe pattern analogous to 
the coherent laser beams and the interference pattern is 
a clear signature of the long range spatial coherence of 
the condensates. The interference pattern generated by 
virtue of the collision of two bright solitons SOI is anal¬ 
ogous to the interference pattern obtained earlier exper¬ 
imentally by Andrews et al. [41 wherein the two con¬ 
densates were separated with a sheet of green light and 


FIG. 13: (a) The time evolution of the trap A (t) = 0.3 —0.09t 2 
corresponding to case (ii), (b) Two soliton interaction form¬ 
ing the interference pattern corresponding to case (ii) for the 
choice a(t) = aoe^°' 15t ^ with aio=0.01, a 2 o=0.8, #io=0.06, 
/3 2 o=0.012, a o =0.02, 0i =0.2, 0 2 =O.l, <h=0.1, S 2 = 0.2, (c) 
Contour plot showing the interference pattern in the confin¬ 
ing regime. 



FIG. 14: (a) The phase difference between the condensates 
corresponding to case (i) and (b) case (ii). 


overlapped in ballistic expansion (switching off the trap) 
while we selectively tune the frequency of the trap A(£) 
in accordance with the scattering length a(t) consistent 
with equation (75). These interference patterns are com¬ 
pletely different from the interference of two independent 
condensates originating from two different traps with or 
without phase [421444] . 

It should be mentioned that the optical traps have 
opened up the possibility of realizing different types of 
temporal variation of the trap frequency A (t) while the 
scattering length can be controlled both by Feshbach res¬ 
onance as well as through the trap frequency A (t). The 
phase change evolution shown in Figs. [l4](a) and[l4| (b) 
gives a measure of the coherence of the condensates. It 
must be added that though the concept of coherence of 
matter waves was already exploited to create atom lasers, 
we do reconfirm the coherent nature of BECs in the intra 
trap collision of bright solitons. 


D. Dynamics of BECs with two and three body 
interactions 
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(a) (b) (c) 



FIG. 15: Density of condensates in the modified GP equation 
with both attractive two-body and attractive three-body in¬ 
teractions for the parametric choice A=0.02, <ao=0.1, /3o=0.9, 
Jx=—1.5, 0i=2.5 for (a) d o =0.4, r = 0, (b)a o =0.4, r = 0.04, 
(c)do=0.45, r = 0.04. 


Relql (a) t = 1 Re|q| (b) t = 1 



FIG. 16: Real part of the order parameter of (a) cubic GP 
equation and (b) modified GP equation with attractive three- 
body interactions. 


(a) (b) 



FIG. 17: (a) Compression and splitting of \q x | of cubic GP 
equation (b) Compression and suppression of splitting of \q x \ 
in the modified GP equation. 


(a) (b) (c) 



FIG. 18: Density of condensates in the modified GP equation 
with attractive two-body interactions and repulsive three- 
body interactions with the same parametric choice as in fig |15| 
for (a) do=0.6, t — 0.04, (b)do=0.8, r = 0.04, (c)do=0.8, 
r = 0.09. 


In equation (50), when the scattering length exponen¬ 
tially increases with time via Feshbach resoncance, the 
density of the condensates also increases [22] [24]. Natu¬ 
rally, a question arises as to how far one can increase the 
scattering length to produce high density condensates. 
Since collapse of the condensates sets in once the density 
exceeds a critical value and one does not expect the col¬ 
lapse of BECs for a true one-dimensional system, there 
is a constraint on increasing the density of the conden¬ 
sates to sustain an effective true one-dimensional BEC 
[22]. This implies that one has to investigate the dynam¬ 
ics of BECs within a safe range of parameters. Hence, 
one has to look for an alternative to generate high density 
condensates retaining the one-dimensionality of the sys¬ 
tem without restricting to a parametric domain. It was 
observed that for a large number of bosons, the repul¬ 
sive three-body interactions can overcome the two-body 
attractive interactions thereby enhancing the stability of 
the condensates m- 

From the above, we understand that two-body interac¬ 
tions alone are insufficient and they have to be suitably 
combined with three-body interactions to extend the re¬ 
gion of stability while increasing the density of the con¬ 
densates in a quasi one-dimensional regime. 

To introduce the new integrable model describing the 
impact of both two- and three-body interactions on the 
condensates, we now consider an additional phase im¬ 
print on the order parameter i/j(x,t) to generate a new 
order parameter q(x,t) as 


in accordance with the following equations; 

e x = -v^H 2 , (so) 

0t = iVr(ipi)* - tp*ip x ) + 4r|t/)| 4 , (81) 

so that the transformed order parameter q(x,t) obeys an 
evolution equation 

x^ 

iQt T Qxx +2a{t)\q\ 2 q+—^-q+4T\q\ A q+4i v /r(\q\ 2 ) x q = 0. 

(82) 

In equation ([82]), a(t)(= doexp(At)) represents attrac¬ 
tive (do > 0) two-body interactions while the real and 
arbitrary parameter r corresponds to the strength of 
three-body interactions assuming that the contribution 
of three-atom collisions to the loss is negligible (by kick¬ 
ing out atoms from the condensate into thermal cloud) 
m-m- This type of engineering the phase imprint to 
generate a new integrable model describing the impact of 
both two- and three-body interactions on the condensates 
is reminiscent of generating solitons by phase engineering 
of BECs of sodium and rubidium m- 


taxi 

2 
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oi jy 
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lq x | 
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q(x, t ) 1 = vp(x, t)e 2l0 ( x,t ), (79) FIG. 19: Asymptotic forms of the density profiles \q x \ of the 

two-soliton solution of (a) cubic GP equation and (b) modified 
where 6(x, t ) is the phase-imprint on the old order param- GP equation (attractive), 
eter ip(x : t). We now engineer the phase imprint Q(x,t) 
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Making use of the gauge transformation approach, we 
obtain the bright one soliton solution of the modified GP 
equation (82) [52] 


q^ = 


A)d 2 4 +2 dsech(xi)exp(i£i), (83) 


Xi 

(X i 


= 2pix — 8 J(a 1 /3 1 )dt + 26 1 , 

= 2a\x — A J (af — fi\ )dt — 2<j>i, 

= a 0 e At , Pi = P 0 e xt , 


(84) 

(85) 

( 86 ) 


and 


0 = — — \/T/3otanh 
a 0 


e At (xA — 2e Xt ao)/3o 
A 


+ *i 


(87) 

where <ao, A), 5i and <A are arbitrary real constants. It 
should be mentioned that the phase imprint 0 given by 
equation (87) is related to the density of the old order 

d80b 


parameter ip (by virtue of equations ( [80] ) and ( |8l[ ) and 
hence the bright solitons given by equation ([83]) are en¬ 
dowed with phase dependent amplitude. It should also 
be added that equation (82) admits bright solitons only 


for attractive two-body interactions (ao > 0) while the 
three-body interactions can be either attractive (r > 0) 
or repulsive (r < 0). 

I. Condensates with attractive three-body in¬ 
teractions (r > 0) 


From the equations (83)-(87), one observes that the 


bright solitons of the integrable modified GP equation 
(82) acquires a kink like additional nontrivial phase rep- 

in comparison with that of cu- 
Figure 15 (a) portrays the 


resented by equation ( |87 
bic GP solitons (r = 0). 
density evolution of the condensates in the absence of 
three-body interactions. When one considers the attrac¬ 
tive three-body interactions in addition to the attractive 
two-body interactions with the strength of both being 
equal, one does not observe any perceptible change in the 
matter wave density as shown in Fig. [l5](b). When the 
strength of attractive two-body interactions is increased, 
one observes a decrease in the density of condensates as 
shown in Fig. [l5] (c). From this, one understands that 
an instability sets in the condensates leading to the ejec¬ 
tion of the atoms resulting in the decrease of matter wave 
density. Hence, to ensure the stability of the condensates 
over a longer interval of time, the strength of the attrac¬ 
tive two-body interactions should be minimum. 

It is also obvious that this additional phase arising 
in the dynamics of the condensates of the modified GP 
equation with attractive three-body interactions evolves 
in space and time as indicated by equation (87) and hence 


one understands that its effect will be more pronounced 
either in the real (or imaginary) part or in the derivatives 
of the order parameter. Figure [l6| (b) shows the effect of 
the nontrivial phase on the matter wave solitons of the 


modified GP equation in comparison with that of cubic 
GP equation shown in Fig. [l6] (a). Thus, it is evident 
that the additional nontrivial phase contributes to the 
compression and rarefaction of the matter wave solitons 
of the modified GP equation. 
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From the spatial derivative \q x \ plotted in Figs. 

(a) and [TT] (b), one infers that as the atoms start accu¬ 
mulating in the lowest quantum state, the matter wave 
density increases leading to a compression of the matter 
wave of the cubic GP equation. In addition to the com¬ 
pression, there is a splitting of the matter wave in the 
cubic GP equation (see Fig. 17 (a)), while this splitting 


of the matter wave has been completely suppressed by 
the attractive three-body interactions as shown in Fig. 

Hoo- 

II. Condensates with repulsive three-body in¬ 
teractions (t < 0) 

For repulsive three-body interactions, the bright soli¬ 
ton of integrable modified GP equation becomes 


qV = 


/3osech(xi)exp(i£i)e( 2 4 +e ) , 


where 0 = 4-yJrfi otanh 2 ( 
a 0 |_ \ 


e t (xA—2e t ao)/3o 




( 88 ) 


with 


f being a real parameter. 

Figure 18 (a) displays the density evolution of the 
condensates for attractive two-body interactions and re¬ 
pulsive three-body interactions. Thus, the addition of a 
small repulsive three-body force contributes to an expo¬ 
nential increase in the matter wave density (as the am¬ 
plitude of the condensates becomes ^|=/3oexp(At/2 + 9) 
from equation (88) in comparison to the density evolu¬ 


tion of the condensates in the absence of three-body in¬ 
teractions shown in Fig fl5|(a). The fact that the con¬ 
densates can hold enormous number of atoms together 
with a small repulsive force means that one can extend 
the region of stability of the condensates in a quasi one¬ 
dimensional regime. When one increases the strength 
of attractive two-body interactions, the density of the 
condensates decreases as shown in Fig. [l8] (b) as com¬ 
pared to Fig. [l8] (a), thereby setting in an instability in 
the condensates. However, this instability can be over¬ 
come by increasing the strength of repulsive three-body 
interactions as shown in Fig. [l8](c), wherein the matter 
wave density again increases enormously in comparison 
with Fig. 18 (b) and attains the magnitude shown in 


Fig. [18] (a). Thus, we observe that the addition of repul¬ 
sive three-body interactions can allow the number of con¬ 
densed atoms to increase enormously and this can hap¬ 
pen even when the strength of the repulsive three-body 
interactions is very small compared with the strength of 
two-body interactions [45]. It should be mentioned that 
eventhough the bright solitons of the modified GP equa¬ 
tion with both two- and three-body interactions (both 
repulsive and attractive) in the expulsive potential are 
unstable again, the extent of instability has reduced com¬ 
pared to that of a cubic GP equation with two-body in¬ 
teractions alone. 
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V. DYNAMICS OF VECTOR BOSE-EINSTEIN 
CONDENSATES 


We know that the behaviour of single component (or 
scalar) BECs is influenced by the external trapping po¬ 
tential and interatomic interaction. Experimental real¬ 
ization of BECs in which two (or more) internal states 
or different atoms can be populated has given a fillip to 
the investigation of multicomponent BECs j53h55] . In 
contrast to the single component BECs, the behaviour 
of multicomponent condensates is much richer because 
of both intra species interaction and inter species in¬ 
teraction. This extra freedom arising by virtue of the 
interaction among the internal states or different atoms 
offers multicomponent BECs several interesting and com¬ 
plicated properties which are not witnessed in single com¬ 
ponent condensates. Far from being a trivial extension 
of the single component BECs, multicomponent conden¬ 
sates exhibit novel and rich phenomena such as soliton 
trains, soliton pairs, multidomain walls, spin-switching 
[561157] and multimode collective excitations. 

Motivated by the above considerations, we investigate 
the dynamics of a two-component BEC in a time de¬ 
pendent harmonic trap described by a two-coupled GP 
equation and deduce the integrability condition for the 
existence of vector bright solitons. We then bring out the 
fascinating collision of vector bright solitons demonstrat¬ 
ing the switching of energy underscoring the longevity of 
vector BECs. We then employ Feshbach resonance man¬ 
agement to manoeuvre the scattering length to further 
enhance the lifetime of two component BECs. 

A. Mathematical Model 

Considering a two-component BEC, the behaviour of 
the condensates that are prepared in two hyperfine states 
can be described at sufficiently low temperatures by the 
two coupled GP equation of the following form [531155] 

in 0 ^ = (-^ y2 + u iM 2 + Uu\H 2 + v) <M89) 
ih^r = (-^ y2 + u 2i\H 2 + U 22 \H 2 + ^ MW) 


be written as 

i ^tt = + 6 nl^i| 2 + & 12|^2| 2 + y X 2 ^ V>1,(91) 

'it - +m« 2 +i* 2 ) « 92 > 

where units for length and time become \ —-— and 

27t/cji_l and ^1,2 is normalized such that f \ipi \ 2 dx = 
1 and j\^ 2 \ 2 dx = N 2 /N 1 . Other parameters are 
defined as: bn = 2an7Vi, 612 = 2miai2W/[(l + 
^2±/^i±)m\, &21 = 2 mia 2 iW/[(l + w 2 ±/uji±)m\, b 22 = 
2a 22 fciVia;2±M±, Ai = u 1x /uju_, X 2 = u 2x /u\± and 
k = mi/m 2 . 

Now, we consider k = 1 (i.e., m 1 = m 2 ) and u)\ = uj 2 = 
uj ( i.e ., Ai = A2) and allow the scattering lengths a^- and 
the strength of trapping potential A 2 (A = cj x /uj±) to vary 
with time. Then, the above equation with t = t/ 2 takes 
the following form (after omitting tilde) 


#it + tlxx + 2(&nO)|'0i| 2 + 6i 2 (i)|V’2| 2 )V’i + X(t) 2 x 2 ipi = 0(93) 
# 2 t + 1p2xx + 2(621 (i) |^ 1 1 2 + b 22 {t)\lp2 1 2 ) 1p2 + X(t) 2 x 2 1 p 2 = Q(04) 

In the above equation, bij(i,j = 1,2) represents the at¬ 
tractive interaction strength and the trap frequency could 
be both confining (A (t) 2 > 0) and expulsive (A (t) 2 < 0). 

When the longitudinal trapping potential is neglected 
(A(t)=0) and bn = b\ 2 = b 2 \ = b 22 = a constant, the 
system becomes the celebrated Manakov model (see Refs. 

[58] to [64]) admitting shape changing collision of vector 
solitons. 

When both intraspecies interaction and interspecies in¬ 
teraction are all equal and time dependent (i.e., bn(t) = 
b 12 (f) = b 2 i(t) = b 22 (t) = a(t)), the above equations ( |93| ) 
and ([94]) take the following form 


iipit + i’lxx + 2<a(t) (|"0i 1 2 + I^2 | 2 )^i + A(i) 2 x 2 ^i = 0(95) 

itp2t + Ip2xx + 2a(6)(|'i/> 1 | 2 + \ip 2 \ 2 )ip 2 + \(t) 2 x 2 ip 2 = 0(06) 

The above coupled GP equation has been investigated 
and shown to admit integrability for symmetric interac¬ 
tion strengths [641166] . 


where the condensate wave functions are normalized 
through the particle numbers TV* = J|^| 2 d 3 r. Uu = 
47r h 2 aa/m and Uij = 2irh 2 aij/m represent intraspecies 
and interspecies interaction strengths respectively with 
dij being the corresponding scattering lengths and m is 
the reduced mass. The trapping potentials are assumed 
to be Vi = rrii[uj 2 x x 2 + uj 2 a (y 2 + z 2 )\/ 2. Further assum¬ 
ing that ujix ojix such that the transverse motions 
of the condensates are frozen to the ground state of the 
transverse harmonic trapping potential, the system be¬ 
comes quasi one dimensional in nature. Integrating out 
the transverse coordinates, the resulting equations for the 
axial wave functions Vq,2 (x,t) in dimensionless form can 


B. Lax-pair and integrability condition 


Equations (|95| and (96) admit the following linear 
eigenvalue problem 


+ £/4> = 0, 


(97) 

= o, 


(98) 

where 4> = (</>i,</> 2,03 ) t and 

( K(t) Qi 

Q2 \ 


U = Qt 

0 

(99) 

\ Q*2 0 - 

-m ) 
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f V n V 12 V 13 \ 
V = I V 2 i V 22 V 23 I , 

V v 31 f 32 ^ 33 ) 


( 100 ) 


with 

vu = -iCW 2 + *r(tK(t) + |QiQt + |Q 2 Q2 
V12 = (r(t)a; — C(t))Qi + -Qix 
V 13 = (r (t)x - C(t))Q 2 + l -Q2x 

V21 = -(mx - amt + l -Q\ x 

V 22 = i((t) 2 -ir(t)x((t) - l -QiQ\ 

V23 = -yQlQl 

v 3 i = -(r (t)x - ami + \Qlx 

V32 = — ~jQ\Q*2 

V 33 = Kit) 2 -iT(t)xC(t) - ^Q 2 Q 2 


Q 1 — 
Q ‘2 = 


1 

7W) 

1 

7W) 


ip 1 a,t)e ir{t)x ‘ 

v> 2 (x,ty T{t)x 


It is obvious that the compatibility condition 
(^ > a;)t = (^ ) t)a: leads to the zero curvature equation 
Ut — V x + [U, V] =0 which yields the integrable cou¬ 
pled GP equation (95) and (96) provided the spectral 
parameter £(£) obeys the following nonisospectral 
condition 



FIG. 20: Switching of matter wave (bright) solitons in two 
component BECs for a(£)=0.01 and A (t)= 0.01+0.0001 t 2 


(96) employing gauge transformation method have the 
following form 

V’i 1) = y^E £ ( 1 )2^ 1 ( i ) sec / l (5» 1 ) e *(-«(104) 


= j £ 2 1)2 AW secft ( #1 ) e ' ( «i+r(*) x 2 ) } (105) 


_(i). 


where 


0r = 4y a 1 p 1 dt + 2xp 1 -2S lt (106) 

Ci = 2 J{a\ - fii)dt + 2 xai - 2 X i, (107) 

with ai = aioe^ r ^ dt \ = /3ioe^ r ^ dt ^ while Si and 
Xi are arbitrary parameters. Gauge transformation ap¬ 
proach can also be extended to generate multi soliton 
solution m • Figure |20| shows that it should be possible 
to transfer (switch) energy from one mode to the other 
in a vector BEC thereby one can enhance the longevity 
(or lifetime) of the bright solitons (or the condensates). 
This can happen irrespective of whether the trap is time 
dependent or independent. 


at) =/J,e-U r ^ dt) (101) 

where fi is a hidden complex constant and T(t) is an arbi¬ 
trary function of time and is related to the trap frequency 


a 2 (t) = r 2 (t) -r'(t). (102) 

Further, equation ( |102| ) which represents the trap fre¬ 
quency A (t) can be related to the scattering length a(t) 
through the integrability condition 

— + 2 a'(t) 2 — A 2 {t)a 2 (t) = 0. (103) 


It should be mentioned that the coupled GP equation 
represented by (95) and (96) is completely integrable for 
suitable choices of the trap frequency \{ t) an d scatter¬ 
ing length a(t) consistent with equation (103). Fo r the 
constant trapping frequency, A (t) = ci, equation ( 103 ) 
yields a(t) = e Clt . The bright solitons of equations (95)- 


C. Feshbach Resonance Management of vector 
bright solitons [68^ 


It is obvious from the investigation of vector BECs that 
the lifetime of the condensates gets enhanced by virtue 
of the switching of energy between the two components. 
It should be mentioned that Feshbach resonance man¬ 
agement can also be employed to drive home the fact 
that vector BECs in a time dependent harmonic trap are 
longlived compared to the condensates dwelling in a time 
independent trap. 

To start with, we now switch off the time dependence 
of the harmonic trap and keep track of the evolution of 
the condensates by manipulating the scattering length of 
the form a(t) = 0.5e~ Q ,25t (show n in Fig. [2l| (f)) con¬ 
sistent with equations ( 102| ) and ( 103| ), thereby render¬ 
ing (r(£) = —0.25) the trap expulsive (shown in Fig. 
(e)). The corresponding density profile of the con- 


21 


densates is shown in the upper panel of Fig. ETl (Figs. 
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|2l| (a) and [2l] (b)). The associated numerically simu¬ 
lated density profile of the condensates employing real 
time propagation of Split Step Crank Nicolson method 
is shown in Figs. [21] (c) and 


21 


(d). From Figs. 


21 


(a)-|2l| (d), we observe that there is a perfect agreement 
between analytical and numerical results. When we dou¬ 
ble the trap frequency (T(t) = —0.5) keeping the trap 
expulsive again (shown in Fig. [22] (e)) and manipulate 
the scattering length through Feshbach resonance of the 
form a(t) = 0.5e _0,5t (shown in Fig. 22 (f)) consistent 


with equations(102) and ( |103| ), the compression (analyt¬ 
ical) the condensates sets in the two modes as shown in 
Figs. [22] (a)-(b). This is again confirmed by the numer¬ 
ical simulation shown in Figs. 22 (c) and [22] (d). When 
we further enhance the trap strength (3.6 times the orig¬ 
inal strength) keeping the trap expulsive again (shown 
in Fig. [23] (e)) and further manipulate scattering length 
a(t) of the form a(t) = 0.5e _0,9t (shown in Fig. [23^f)), 
one observes the onset of collapse of the condensates as 
shown in Figs. [23] (a)-(b). Again, the numerically simu¬ 
lated condensates shown in Figs. 23 (c) and|23](d) match 
with Figs. 23 (a) and 23 (b) respectively. 

To enhance the stability of the condensates, we now 
switch on the time dependence in the harmonic trap keep¬ 
ing it expulsive (shown in Fig. [24] (e)) of the form T(t) = 
—0.25£ and manipulate the scattering length through Fes¬ 
hbach resonance of the form a(t) = Q.5e ~°‘ 12 5t sho wn i n 
Fig. [24] (f) consistent with equations (102) and (103), 
the corresponding density profile is shown in Figs|24|(a)- 
(b). This is confirmed by the numerical simulation of the 
condensates shown in Figs. [24] (c) and 24 (d). When we 
make a 20 fold increase in the expulsive trap frequency 
(shown in Fig. [25] (e)) and accordingly employ Feshbach 
resonance to choose a(t) = 0.5e~ 2-5t (shown in Fig. 25 


(f)) consistent with equations ( 102| ) and ( |103| ), the cor¬ 
responding density profile shown in Figs. |25[a)-(b) ex¬ 
actly match with the numerically simulated condensates 
in Figs. [25] (c) and [25] (d). 

When we further enhance the time dependent expul¬ 
sive trap (shown in Fig. [26] (e)) frequency T(t) by 100 
times the original strength and manipulate scattering 
length (shown in Fig. [26] (f)) through Feshbach reso¬ 
nance (in accordance with equations ( |102| and |103| )) , one 
never sees an abrupt increase in the density of the con¬ 
densates as shown in Figs. [26] (a)-(b). Thus, we observe 
that any further increase of the expulsive trap frequency 
T (t) and manipulation of scattering length through Fesh¬ 
bach resonance does not significantly increase the density 
of the condensates eventhough the attractive interaction 
strength increases rapidly. In other words, the conden¬ 
sates in the time dependent expulsive harmonic trap con¬ 
tinue to remain stable for a reasonably large interval of 
time even for large attractive interactions. Again, an¬ 
alytical results synchronize with numerical simulations 
shown in Figs. [26^c)-(d). 

Thus, we observe that the vector BECs in a time de¬ 
pendent expulsive trap are more long lived and the life¬ 
time can be enhanced by Feshbach resonance manage- 


M 


M 



FIG. 21: Upper Panel (a)-(b): Density (Analytical) profile 
of the condensates (bright solitons) in the time independent 
expulsive trap for T(t)= -0.25 and a(t) = 0.5e _O 25t ; Mid¬ 
dle Panel (c)-(d): Numerically simulated density profile for 
r (t)= -0.25 and a(t) = 0.5e _o ' 25t ; Lower panel(e)-(f): Trap 
strength and binary interaction 


ment compared to the condensates in a time indepen¬ 
dent expulsive harmonic trap. We also wish to point out 
that we are able to stabilize the condensates in an expul¬ 
sive time dependent harmonic trap for attractive interac¬ 
tions, where the condensates usually get compressed and 
collapse subsequently. It is also pretty obvious that the 
scalar counterpart of the condensates stabilized through 
Feshbach resonance (shown in Figs. 24]|26 ) will collapse 
immediately during time evolution. 


VI. STABILIZATION OF BRIGHT SOLITONS 
IN WEAKLY COUPLED BECS 

A. Impact of weak time dependent Rabi 
coupling 

We also emphasize that the enhancement in the life¬ 
time of the two component BECs arises by virtue of in¬ 
traspecies and interspecies interaction. It should be men¬ 
tioned that one can make two component BECs more 
longlived by the addition of weak time dependent or space 
dependent coupling forces. 

We now consider a spinor BEC composed of two hy- 
perfine states, say of the |F = 1, m/ = —1 > and 
|F = l,ra/ = l> states of 87 Rb atoms [69] confined 
at different vertical positions by parabolic traps and cou¬ 
pled by a time dependent coupling field. 
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FIG. 22: Upper panel (a)-(b): Compression (Analytical) 
of the condensates in the time independent expulsive trap 
for T (t)= -0.5 and a(t) = 0.5e _o ' 5t ; Middle Panel (c)- 
(d): Numerically simulated density profile for T(t)= -0.5 and 
a(t) = 0.5e _o ' 5t showing the compression of BECs; Lower 
panel (e)-(f): Trap strength and binary interaction. 


FIG. 24: Upper panel (a)-(b): Density (Analytical) profile 
of the condensates by switching ON the time dependence in 
the expulsive trap for T(t)=-0.25 t and a(t) = 0.5e -o ' 125t ; 
Middle Panel (c)-(d): Numerically simulated density pro¬ 
file for r(t)=-0.25 t and a(t ) = 0.5e _o ' 125t ; Lower panel 
(e)-(f): Transient trap and Interaction strength 



FIG. 23: Upper panel (a)-(b): Onset of collapse (Analyt¬ 
ical) of the condensates for T(t)= -0.9 and a(t) = 0.5e _O-9t 
in an expulsive time independent trap; Middle Panel (c)- 

(d) : Numerically simulated density profile for T(t)= -0.9 and 
a(t) = 0.5e _o ‘ 9t showing the collapse of BECs; Lower panel 

(e) -(f): Trap strength and binary interaction 


FIG. 25: Upper panel (a)-(b):Density (Analytical) pro¬ 
file of the condensates in the time dependent expulsive trap 
for r(t)=-5 t and a(t)=0.5e -2 ' 5t ; Middle Panel (c)- 
(d): Numerically simulated density profile for T(t )=-5 t and 
a(t)=0.5e _2 ‘ 5t ; Lower panel (e)-(f): Transient trap and 
Interaction strength 




































18 



FIG. 26: Upper panel (a)-(b): Density (Analytical) pro¬ 
file of the condensates in the time dependent expulsive trap 
for r(t)=-25 t and a(t) — 0.5e _12 ' 5t ; Middle Panel (c)- 
(d): Numerically simulated density profile for T(t)=-25 t and 
a(t ) = 0.5e _12 ' 5t ; Lower panel (e)-(f): Transient trap and 
Interaction strength 



FIG. 27: Density of the condensates for a(t) = 0.5t, = 0.3 
and T (t) = O.lt X 10“ 2 without coupling. 


We assume the condensate to be quasi-one-dimensional 
(cigar-shaped). Then, in the meanfield approximation, 
the system is described by the coupled GP equation Uoj 


iipit + ipixx + 2a(i)(|^i| 2 + |V , 2| 2 )V’i + \{t) 2 x 2 -^i 

+ = 0, (108) 

ifp2t + 1p2xx + 2a(i)(|V’l | 2 + \'lp2\ 2 )lp2 + \(t) 2 X 2 1p2 

+ iG{t)^2 + ^(£)Vh = 0? (109) 


In the above equation, v(t) denotes the coupling be¬ 
tween the two condensates while G(t) accounts for the 
feeding (loss/gain) of the condensates from the thermal 
cloud. Infact, the impact of adding a linear coupling 
to the Manakov model was already demonstrated m 
and the resulting dynamical system was sho wn be inte¬ 
grate. The linearly coupled GP equations (109) have 



FIG. 28: Impact of weak coupling on the con den sates with 
u(t) = 0.11 with the other parameters as in Fig|27| 



FIG. 29: Stabilization of the condensates by tuning the scat¬ 
tering length for a(t) = 0.5 X 1/2 1 with the other parameters 
as in Fig. [28] 


also been investigated and the concept of domain walls 
m and symmetry breaking of solitons m was explored. 
Equations (109) have also been investigated for v(t) = 0 
[66l 67. and the dynamics of vector BECs has been ana¬ 
lyzed. The explicit forms of one soliton solution of equa¬ 
tions (109) can be written as m 


ip[l) =/3 1 (t)sech(6i)e' l< ' G+r(t) x 2 )^ ( 110 ) 
V a (t) 

•02 1 ) = Pi(t)sech(0 i) e * ( “« 1+r(t) V), (HI) 


where 

6i = 8 J ai(t)/3i(t)dt + 2x/3i(t) — 25i, (112) 

£i — 4 J (a\(t) 2 — /3i(t) 2 )dt + 2xa\(t) — 2xi,(113) 

with ol\ (t) = ai^e~ f 2r P) 2<it , p 1 (t) = Aoe - J' 2 r P) 2rft 
while Si and xi are arbitrary parameters with e 1,2 as 
coupling parameters. Thus, it is obvious from equation 
( |110| ) and ( |111| ) that the amplitude of the bright solitons 
depends on the temporal scattering length a(t) and trap 
frequency T(t). (Pi(t) varies exponentially with T(t)). 
The fact that the trap frequency T(t) varies exponen¬ 
tially with the time dependent coupling coefficient v(t) 
(see equation (9) in Ref. [74]) indicates that the density 
of the bright solitons or the condensates could build up 
in a very shot span of time. This means that the dynami¬ 
cal system could enter into the domain of instability very 
quickly. 

The profile of the density of the condensates without 
the linear Rabi coupling is shown in Figs. [27] (a) and (b). 
When the two condensates are coupled together, one ob- 
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(a) (b) 



FIG. 30: Stabilization of the condensates by finetuning the 
coupling coefficient for v(t) — O.lt/5 with the other parame¬ 
ters as in Fig. [28] 


(a) <b) 



FIG. 33: Retrieval of intramodal inelastic interaction of the 
condensates for v(t) — O.OltXlO -2 with the other parameters 
as in Fig. [32] 



FIG. 31: Intramodal Inelastic collision of the condensates 
without coupling for a(t) = 5, F(t) = O.lt X 10 -2 , = 0.89z 

and =0.6 


serves an abrupt increase in the density of the conden¬ 
sates as shown in Figs. 28 (a) and (b). This suggests 
that coupling the two condensates leads to an instability 
in the dynamical system. However, this instability can 
be overcome either by changing the temporal scattering 
length through Feshbach Resonance as shown in Figs. [29] 

(a) and (b) or by finetuning the time dependent coupling 
co-efficient as shown in Figs. 30 (a) and (b). 

The gauge transformation approach m can be 
easily extended to generate multisoliton solutions. 
From the collisional dynamics of the condensates, one 
observes an intramodal inelastic interaction of bright 
solitons in the absence of weak coupling as shown in 
Figs. j3lja) and (b). The addition of a weak time 
dependent coupling disrupts the intramodal inelastic 
collision of the condensates as shown in Figs. [32^ a) and 

(b) . However, one can suitably finetune the coupling 
coefficient v(t) to retrieve the intramodal inelastic col¬ 
lision of the condensates as shown in Figs. [33] (a) and (b). 


(b) 



FIG. 32: Impact of coupling on the collision of condensates 


for v(t) — O.Olt with the other parameters as in Fig. 31 


B. Influence of weak spatially dependent 
coupling 

To understand the effect of weak spatially depen¬ 
dent coupling, we now consider cigar-shaped (quasi-one- 
dimensional) BECs composed of two hyperfine states of 
87 Rb atoms m confined by the parabolic trapping po¬ 
tential subject to a spatially dependent force v(x). The 
system is now described by the following coupled Gross- 
Pitaevskii (GP) equation [70 within the framework of 
the mean-field description as 

iipit + ipixx + 2a(t)(\4>i\ 2 + \4> 2 \ 2 )4 >i+ 

+ iG(x,t)ipi + v(x)ip 2 = 0, (114) 


ilp2t + ^2x 


2 a (^)(IV’l| 2 + IV’2n)V’2 + \(t) 2 X 2 Tp2 


+ iG(x,t)ip 2 + H x )^i = 0 > 


(115) 


In the above equations, i/(x) denotes the weak spatial 
coupling between the two components while G(x,t) ac¬ 
counts for the feeding which can be carried out by means 
of a reservoir filled with a large amount of the conden¬ 
sate to which the trap is weakly coupled [76] . In fact, 
effects generated by adding the linear coupling to the 
Manakov model were already studied [? ] and the re¬ 
sulting sy stem was found to be integrable. Equations 
(114) and (115) have been investigated for v = 0 


and the ensuing dynamics of the vectorial BECs has been 
analyzed. When the trapping potential and (weak) spa¬ 
tial coupling between the two condensates are neglected 
(A (t) = v{x) = 0) and the scattering length a(t) becomes 
a constant (neglecting gain/loss G(x,t)) then, the dy¬ 
namical system becomes the celebrated Manakov’s model 

m 

The exp licit f orms of bright soliton solution of equa¬ 
tions (114) and (115) can be written as m 


vP = 




ePft (i)sech(6> 1 )e i (-€ 1+r W / , (116) 


V>P = - 7 F^4 1) AWsech(0 1 )e i( - ?1+r(t U^), (117) 

v a (c 
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FIG. 34: The density of the condensates for a(t) — 0.5£, — 

0.3 and T(t) = f 0.005t dt in the absence of the coupling. 




FIG. 35: The impact of weak linear coupling u(x) = 0.5x on 
the condensates with the other parameters as in Fig. |34| 


where 


6»! = 2 J Pidx + J(8a 1 p 1 -4a 1 )dt-25 1 , 

6 = -2 / a\dx — 4 J (c*i + 4 iff - 4 Pi)dt 


(118) 

nm 


FIG. 36: The evolution of bright solitons (given by equation 
(116)) in the presence and absence of the linear coupling in a 
confining trap for T(t) = f 0.005 1 dt. 



with aq(t) = aqo^ -27 = /3ioe -27 W, where 

and xi are arbitrary parameters and e\,2 are coupling 
parameters. 

Thus, it is obvious from equations ( |116| ) and ( |117[ ) that 
the amplitude of the bright solitons depends on the time- 
modulated scattering length a(t) and /3i(t), which implic¬ 
itly depends on v{x) (by virtue of equation (7) in Ref. 
M)- The density profile of the condensates shown in 
Fig. [34] (a), (b) and Fig. [35] (a), (b) and the time evolu¬ 
tion of the bright solitons shown in Fig. [36] demonstrate 
that, in the absence of spatial coupling, the density of 
the condensates grows with time while they remain lo¬ 
calized around the same point. The introduction of the 
linear coupling not only stretches the wave packet, but 
also shifts the center of the localization of the wave packet 
as shown in Fig. [36] 

The explicit forms of the dark-soliton solution can be 
written as 


FIG. 37: The evolution of a dark soliton (given by equation 
(120)) in the presence and absence of the spatially-dependent 
linear coupling in a confining trap for T(t) = f 0.005 1 dt. 



FIG. 38: Density profiles of bright-bright soliton interaction 
for the parametric choice a(t) = — cos(v / 2t), T(t) = f O.OOlt 
dt, is(x) = [log(—cos (2a:))] 2 , — 0.89z, and — 0.6 in 

an expulsive trap. 


V4 1} = -^=£ i i ) Pi(t)tanh{9 1 )e i ^ 1+r ^^J' xdx ^ (120) 
V4 1} = Pi(t)tanh{0i)e i ^ 1+r< ' x ’ t ^ xdx ) (121) 

\A<t) 

The time evolution of the dark solitons indicates that one 
observes a similar impact of the linear coupling (a shift 
of the center of the localization as shown in Fig. [37] and 
stretching of the matter wave packet) in the dark solitons 
as well. 



FIG. 39: Contour plots for the bright-bright soliton interac¬ 
tion shown in Fig. [38] 
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FIG. 40: The interference pattern arising from the dark-dark 
soliton interaction. 


(a) (b) 



40 40 


FIG. 41: Density profiles for the dark-bright soliton in¬ 
teraction in a transient trap for the parametric choice 
v(x) — [tan(0.03x)] 2 + 0.0075, a{t) — O.OOlt, T(t) = 
— f 0.05 te~°' 005t dt, = 0.89z, and — 0.6. 


FIG. 42: The transient trap for T(t) = — f 0.05 te °' 005t dt. 


(a) (b) 



FIG. 43: Contour plot of the the dark-bright soliton interac¬ 
tion shown in Fig. [4T| 


Another interesting observation is that the addition 
of the spatially-dependent linear coupling does not con¬ 
tribute to the growth of the condensates unlike the time- 
dependent coupling [74]. Instead, the spatial coupling 
merely stretches the wave packet. This means that one 
can easily stabilize the vectorial BECs by the introduc¬ 
tion of the spatially-dependent linear coupling. One can 
easily extend gauge transformation approach to generate 
multiple dark solitons. 

The shift of the center of the matter-wave packet, 
which is induced by the spatially-dependent coupling, 
can be suitably exploited to generate matter-wave inter¬ 
ference pattern which was addressed theoretically m ■ 
The density profile for bright-bright soliton interaction is 
shown in Figs. [38] (a) and (b). From the density pro¬ 
file, one observes that the matter wave density reaches 
the maximum value periodically and the central fringe is 
always bright in every interference pattern. This means 
that the bright - bright solitons constructively interfere 
at the centre giving rise to maximum intensity. One ob¬ 
serves alternate bright and dark fringes in the pattern. 
The contour plot shown in Fig. [39] confirms this obser¬ 
vation. 

The contour plot corresponding to the density profile 
of dark-dark soliton interaction is shown in Fig. [40] In¬ 
terference pattern shown in Fig. [40] indicates that the 
central fringe is always dark indicating that the dark- 
dark solitons destructively interfere at the centre. Again, 
one notices alternate dark and bright fringes. 

Choosing a transient trap of the form shown in Fig. [42] 
the density profile of dark-bright interaction is shown in 
Fig. [4l] From the figures, we understand that when the 


solitons make a transition from the expulsive region to 
the confining region, the dark and bright solitons inter¬ 
fere with each other in the confining domain. When the 
trap gets back to the expulsive domain again, the inter¬ 
ference pattern disappears while the density of the con¬ 
densates gets compressed. The above interference pat¬ 
tern is quite identical to the one observed in Ref. [40] . 
From the interference pattern shown in Fig. [43] one ob¬ 
serves that the central fringe is always bright indicating 
the maximum intensity. This occurs because the ampli¬ 
tude of bright solitons predominates over that of the dark 
solitons. It must be mentioned that the spacing between 
the bright (or dark) fringes in the interference pattern 
represented by Figs. [39] and [40] gives a measure of the 
coherence of matter waves and is a characteristic feature 
of the atomic system under investigation. The spacing 
between the bright (or dark) fringes in the interference 
pattern shown in Fig. [43] observed in the confining trap 
has narrowed down as the condensates dwell in the con¬ 
fining trap for a very short interval of time. 


VII. FUTURE DIRECTIONS 

In the present review, we have analytically proved that 
vector BECs are longlived compared to scalar BECs. 
Eventhough this paper has comprehensively reviewed the 
recent developments from the perspective of integrabily 
in scalar and vector BECs in a time independent/ depen¬ 
dent harmonic trap, there are several unexplored territo¬ 
ries in the domain of ultra cold atoms. The identification 
of an integrable model for a two component BEC (either 
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for a BEC comprising of hyperfine states of the same 
atom or two different atoms) characterized by three dis¬ 
tinct scattering lengths, namely two distinct intraspecies 
scattering lengths and an interspecies scattering length 
has eluded our observation so far. Eventhough weak cou¬ 
pling of the condensates generates high density matter 
waves and one can somehow stabilize them, the ques¬ 
tion of nonlinearly (strongly) coupling the condensates 
remains unanswered. These investigations will certainly 
enable us to penetrate deep into the dynamics of ultra 
cold atoms with precision. Investigations along these di¬ 
rections are in progress and the results will be published 


later. 
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